Excited DeSitter brane worlds localized by a kink. 
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We reconsider, in five-dimensional space-time, the issue of thick brane localized in the extra 
dimension by a kink formed by a scalar field. The localization is achieved by a sine-Gordon potential. 
Apart from a fundamental brane (discovered by Koley and Kar [lj]) where the scalar field is a 
monotonic function of the extra dimension), we show that a series of new solutions exist as well, 
labelled by the number of zeros of the scalar field. These solutions are regular, localized on the 
brane and mirror symmetric with respect to the extra dimension. They form a tower of "excited 
branes". The study of some perturbations of the solutions reveals that the new solutions are not 
stable. Finally, fermions are coupled to the scalar field by means of a Yukawa potential and their 
localization in the background of the new solutions is examined. It turns out that the excited branes 
can localize left and right chiral fermion either on the brane and/or in the bulk but close to the 
brane. 
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I. INTRODUCTION 



One of the main problem in the theoretical description of the Universe is the cosmological constant problem 
0. One of the most promising attemp to solve it is based on the " brane- world" models p|; the basic idea 
is that our Universe is represented by a 3+1-dimensional subspace (a three-brane) embedded into an higher 
dimensional space-time (the bulk). In this context, the first challenge is to reconciliate the extra dimensions 
with the fact that low energy physics is very well described by a conventional 3+1-dimensional space-time in 
which numerous laws (e.g. Newton's law) are tested with a very high degree of accuracy. It is nevertheless 
challenging to emphasize the Einstein equations in higher dimensional space-times and to study and classify 
the possibly relevant solutions. One way to attack the various problems consists in trying to localize the 
matter fields and gravity on the brane in the case where the extra dimensions are non compact 0, Q ; another 
way would be to have solutions where the geometry of the extra-dimensions naturally comes out to be a 
compact manifold 

Here we will follow the ideas developped in several papers [6 8] and consider branes which are supported 
by appropriate topological deflects living in the extra-dimensions. In these papers, the brane is supposed 
to be Minskowski space-time (no vacuum energy) or de Sitter space-time with a cosmological constant (or 
brane tension) leading to an inflating brane. No bulk cosmological constant is assumed. In contrast, in 
Q the authors take advantage of a cosmological constant in the bulk : A4 + „ ^ and manage to enforce 
the localisation of gravity by a fine tuning of A4+ n in function of the other constants characterizing the 
topological defects. 

Assuming only one extra dimension, a lot of work has been done about thick brane scenarios [9l4l^| where 
five-dimensional gravity is coupled to one or more scalars and an appropriate potential leading to localized 
solutions in the extra dimension. The relation of gravity localisation with different types of potential for one 
scalar field was analyzed namely in [l3| and in [T3 - [l6j ; here the relevant potential derive from a superpotential 
allowing for explicit solutions to the field equations. The Sine-Gordon soliton potential was used in [l[. In a 
series of papers, two scalar fields were used to construct a so called Bloch brane [l7l - fl9| which is also based 
on a superpotential. 

Here, we will reconsider the existence of branes within the sine-Gordon model. One feature of the brane 
solution found in [l[ is that it exists only for a particular relation between the five dimensional Newton 
constant and the bulk cosmological constant A. In this paper, we reinvestigate the model of [l[ and study 
the field equations numerically for generic value of the constant A. Our numerical result strongly suggest 
that the Koley-Kar brane-world is the first solution of a tower of solutions characterized by the number of 
nodes of the scalar field. All these solutions are mirror symmetric with respect to the four dimensional brane 
and exist for specific values of the cosmological constant (like a spectral value while the five dimensional 
Newton constant is fixed. For the other values of A, the warp factor become singular on one side of the 
bulk space. The excited solutions have important consequences on the localization of masless fermion on the 
brane. In particular, both chiral component of the fermion can be localized on or close to the brane while it 
is know that only the left component is localizable in the background of the Koley-Kar fundamental solution. 
We also checked that the tower-like solutions of the type studied here are also present with other type of 
potential admitting kink in the flat space, e.g. the so called Xcj) 4 potential. The model and the equations 
are presented in Sect. II. In particular we consider the case where the four dimensional brane is an Einstein 
space with a independant cosmological constant. The various type of solutions are discussed in Sect. Ill 
and Sect. IV is devoted to the localization of the fermions. Finalize, we analyze the stability property of the 
first few elements of the tower in Sect. V and give some conclusions in Sect. VI. 



II. THE MODEL AND THE EQUATIONS 



Along with Koley and Kar [l| (and more recently [2(|), we consider the action 

S = J d 5 x^[^(R - 2 A) - ±g MN d M ct>d N ct> - V (</>)] (1) 

where k§ = 8ttG^ with G5 the 5-dimensional Newton constant and A is the 5-dimensional cosmological 
constant. The real scalar field <j> interact through the potential V ((/>). In this note we have considered both 
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the sine-Gordon potential and the quartic double potential 

VU) = p(l + cos( — )) or V(4>) = p($- - l) 2 (2) 

q q 

where p, q are real constants. These potentials lead qualitatively to the same results but we will discuss 
essentially the case of the sine-Gordon potential. We are interested in solutions of the field equations 
associated to the action (JIJ with a metric of the form 

ds 2 =e- 2A M(dsl)+dy 2 (3) 

where the extra coordinate is denoted y and the function exp(— 2A(y)) is the warp factor (note that the 
function A(y) in our notation differs by a sign from the one of 20]). The scalar field is assumed to depend 
on y only : <ft — f(y)- The 4-dimensional space whose metric is denoted ds 2 and is assumed to be an Einstein 
space, i.e. such that 

a® - (4) 

where H 2 is a constant which can be negative. A real (resp. negative, nul) H 2 corresponds to deSitter (resp. 
Anti de Sitter, Minkowski) 4-dimensional space-time but the metric can actually describe any Einstein 
space. 

The relevant components of the 5-dimensional Ricci tensor are then 

R° = R\ = Rl = R\ = 3e 2A H 2 + A" - 4(A') 2 , R\ = 4A" - 4{A') 2 (5) 
from now on, the prime denotes the derivative with respect to y. The field equations read 

A" = ^(/') 2 +ff 2 ex P 2A , k = k 2 (6) 

A' 2 = y 2 {J') 2 - \kV{J) - | +ff 2 ex P 2A (7) 
dV 

f" = AA 'f' + ^j (8) 
of 

These equations depend in principle on five constants : the parameters of the potential p, q the gravity 
parameter k and the bulk and the brane cosmological constants, respectively A, H . The parameter q defines 
the scale of the scalar field and can be set to one. Independently, p can be set to a fixed value by a suitable 
rescaling of y, A and of k. In the following we adopt the normalisation q = 1, p = 1/4. Using (j6|) and dHJ to 
determine A(y),f(y), the cosmological constant A then appears as a "constant of motion" through Eq.([7]). 
As a consequence, only k, H need to be specified to construct a solution. 



III. ANALYSIS OF THE SOLUTIONS 



A. Case H = 



The equations ©,© have to be solved with appropriate boundary conditions. Without loosing generality, 
the coordinate y can be translated in such a way that 0(0) = 0; we want the sine-Gordon field to approach a 
vacuum for y — > oo, so we require (p(oo) — n/2. In addition, the condition A(0) — can be set by a suitable 
normaliation of t , the brane is then located at y = 0. This leaves the parameter A'(0) = C to be chosen in 
order to specify a boundary value problem. 

The hypothesis of a symmetric space time under the reflection y — > —y is fulfilled if C = 0. However, 
assuming for a while space-time limited to the subspace with y > 0, the equations can be integrated for 
y € [0, oo] with arbitrary values of C. The boundary conditions are then 



A(Q)=0 , A'(Q) = C , /(0)=0 , /M = - 



(9) 
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FIG. 1. The functions A, A', f, f for « = 1, p = 0.25 and G = 0.1 



For C = 0, the explicit solution of Koley, Kar is recovered 



-4(2/) = g lncosh(fcy) , f(y) = 2 arctan(exp %) - ^ (10) 



with 



, 2 12p A 



fc^ = — ^— , A = --^ (11) 

4k + 3 4k + 3 y ' 

In particular, the bulk cosmological constant is fixed by the five dimensional coupling constant and the 
parameter p. 

For C > 0, the system of equation above has, to our knowledge, no explicit solution. We therefore solve 
the equation by a numerical method. We used a collocation method for boundary- value ordinary differential 
equations, equipped with an adaptive mesh selection procedure [22j]. Our solutions were constructed with 
a relative error of order 10 -8 . The solution corresponding to k 5 = 1 and C = 0.1 is presented on Fig. Q] 
Continuing the numerical solution on the negative part of the y-axis, the numerical results show that, for 
generic values of C, the solution develops oscillations and then becomes divergent when y is small enough. 
Accordingly, the solution cannot be continued of the full axis y € [—00,00]. The analysis reveals, however, 
that for fine tuned values of the parameter C, the solution has f(y c ) — A'(y c ) — for some y c . For these 
values of C and y c , the solution can be continued into a regular, mirror symmetric solution (after a translation 
by Vc) on the full axis y. 

We believe that solutions with an arbitrary number of oscillations in the core exist, forming an infinite 
tower labelled by N, the number of nodes of f(y) on y £ [0,oo] (corresponding to a total on 2N — 1 
nodes). We constructed the first three elements of this tower; the profiles are presented on Fig. [5] (A, f 
on the left part, A',f on the right part). Several parameters characterizing the brane solutions, namely 
the bulk cosmological constant and the value |/'(0)| are represented as functions of K5 on Fig. [5] for the 
first three s olution s. The curves corresponding to N — 1 is known analytically; the A function is such that 

A' (00) = v/^A/a 

The pattern of solutions of the sine-Gordon equation (i.e. in the probe limit K5 = 0) suggests that the 
new solutions, found for k > 0, are nothing else than the first elements of a tower of gravitating solutions 
labelled by the number of nodes of the scalar function. We now present arguments supporting this statement. 
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FIG. 2. The metric function A and scalar function / for the solutions with one, two and three nodes(left), the 
corresponding derivatives (right) 




FIG. 3. The parameters |/'(0)| and A as functions of K5 for N = 1, 2, 3 



With our conventions (q — l,p = 1/4 and assuming the kink to be centered at the origin: /(0) = 0), the 
sine-Gordon kink has /'(0) = 1, /(oo) = ir/2. For the following, let us denote a "half-kink" the part of this 
solution between y = and y = 00. Integrating the sine-Gordon equation with /(0) = 0, /'(0) = s with s < 1 
leads to oscillating solutions; these are periodic functions whose zeros (which we note y = y n ) are such that 
f'iVn) — ±s- Accordingly, configurations can be constructed by "gluing" a half-kink (suitably translated) 
at one of the zeros of a periodic solution corresponding to a value of s arbitrarily close to s = 1. Such 
configurations are off course not differentiable at y = y n but the coupling to gravity somehow regularizes 
them and leads to the family of regular solutions discussed above. Our numerical results strongly confirms 
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this interpretation. 

Replacing the sine-Gordon potential by the double well quartic potential given in ([2]) leads to the same 
pattern. The occurence of the " quasi-solution" configurations in approaching the probe limit makes, however, 
the numerical analysis rather involved for K5 << 1 for the two potentials. 

B. Case H 2 < 

Let us pose H 2 = —K 2 . In this case, we could not find explicit solutions for generic values of k. For 
«5 = 0, the equation for A can be solved, leading to 



A(y) = ilog(l 



tanh 2 (#?/)) 



exp(-2A) = cosh 2 (iJy) 



(12) 



This function is such that A = —K 2 y 2 /2 for y << 1 and A = — \K\y for y >> 1. Accordingly, the warp 
factor exp(— 2 A) increases for large y and the brane cannot be localized. Solving the equations numerically 
for K5 > 0, reveals that, asymptotically, the warp factor keeps diverging for large y. Several plots of such 
solutions are presented in Fig. 4. 




(a) 



y 




(b) 



y 



FIG. 4. Left: The metric function A and scalar function /, /' for k = and three values of H 2 . Right: the metric 
function A (solid) and the derivative (dashed) for k = 1 and three values of H 2 . 



Case H > 



In this case also, we could not find explicit solutions for generic values of n. For n = 0, the equation for 
A can be solved, leading to 



1 



A(y) = -log(l + W(/fy)) 



exp(-2A) = cos 2 (Hy) 



(13) 



It is such that A(y) = H 2 y 2 /2 + 0(y 3 ) for y << 1. The warp factor vanishes y = tt/ (2H), while the compo- 
nents of the Riemann tensor remain finite. Solving (numerically) the equation for f(y) in this background 
reveals that the scalar field does not stabilize to a constant for y — > ir/(2H). 

Investigating the equations numerically for n > reveals that the peculiar behaviour of the k = limit 
(in particular the vanishing of the warp factor at a finite value extra dimension coordinate y) persists. 
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IV. LOCALIZATION OF THE FERMIONS ON THE BRANE 



The coupling of the Dirac action to the brane is usually realized by means of the fermionic action 

S F = I d 5 x^($T M D M ^ - rmpip - ^F{4>)^ • (14) 



The parameter r\ represents the coupling of the fermion to the scalar held. The function F(<fi) is usually 
chosen as F{<jj) = </> leading the a conventional Yukawa potential, (the alternative choice F((j>) = sm(<j>) is 
also used (ll. |20||). The matrices T M are the Dirac matrices in curved space. To solve the underlying Dirac 
equations, one usually adopts a coordinate z defined according to dz = exp(A(y))dy. This transforms the 
line element into a conformally flat metric 

ds 2 = e-^W^dafda? + dz 2 ) (15) 

With this coordinate, the Dirac matrices are given by T M — (e A ^^ 1 e A ^ 5 ) where 7 A ',7 5 are the usual Dirac 
matrices in flat space. The fermion function is usually decomposed according to 

ij}(x, z) = 2J i>L,nOiLn{z) + X! ^R,n a Rn( z ) ( 16 ) 

n n 

with the chiral spinor ipL,n = — 7 5 V J £,n an d ipR,n = 7 5- 0-R,n and the sum n run over the spectrum of the 
4-dimensional solutions. 

After some standard manipulations (see e.g. in [H,[I(|), the effective potential determining the ^-dependent 
prefactors a.L n (z), cti, n (z) for the chiral components of the Dirac spinors are given by 

%M=e-W-^+^U) ( 17 ) 

VR(») = e-"tf? + 1^-^)1^ (18) 

It is known that the fundamental Koley-Kar brane leads to a vulcano shaped potential Vl which can localize 
the fermion on the brane while Vr is bell shaped, admitting no bound states. Accordingly, the right handed 
fcrmions cannot be localized on the brane. 

Owing on the new excited branes constructed in the previous section, it is natural to study the shape of 
Vl and Vr corresponding to the excited solution. This is illustrated by Fig. [5] It turns out that the shape 
of the potential corresponding to the exited brane is quite different from the one of the fundamental brane 
and presents some interesting features. On Fig|S]the potentials Vl,r corresponding to k = 1, r\ = 1 and the 
N=l and N=2 solutions are superposed. For the first excited brane, it turns out that Vr presents a negative 
valley on the brane (i.e. for z = 0) and at least one bound state should exist there. In contrast, the potential 
Vl presents a local maximum on the brane and a local minimum at some finite value y = ±y m . It means 
that left handed fermions could be bounded in this region; the corresponding fundamental wave function 
should be localized symmetrically in the two valleys. Decreasing the Yukawa coupling constant k, it turns 
out that the potentials Vl and Vr both present a local minimum on the brane; accordingly both species of 
fermions could be localized. 

Taking the next excitation, reveals that the corresponding effective potentials develop more and more local 
minima in the bulk, as shown on Figj6j 

V. STABILITY 
A. Setup 

We consider the following ansatz for the metric and matter field perturbations: 

ds 2 - e-^+^^n^dx" + e-^ t ^dy 2 + eC 1 (t,y)dtdy , f(t,y) =f + efx(t,y), (19) 
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FIG. 6. The effective potentials for the first (left) and second (right) excited branes for r) = 1 and two values of the 
: k parameter: k — 1, 0.003 for the left and k = 1, 0.006 for the right graph. 



where e is a small bookkeeping parameter and A, /o are solutions to the background equations given in the 
previous section (note that, for convenience, the function / of the previous section is here renamed /o). The 
form P^|) is the most general ansatz compatible with the symmetries of the background spacetime. The 
stability issue of thick branes in five dimensions has been adressed in [2l[ and we basically follow the same 
approach. 

Assuming, as usual, the fields to depend on time through a factor e , the nontrivial gravity perturbation 
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equations reduce to 

3A' (8A[ - B[) + B x (\2A a - K ff) + 2 (-3A'( + nfif[ + k/xV'C/o)) = 0, 

-24A'A[ - \2B x A n + Ziuje^dA 1 + 6uj 2 e 2A A 1 + nB x f 2 + 2k/' f[ - 2Kf 1 V'(f ) = 0, 

WiA' + 6A[ - 2Kfifo = °- (2°) 

The perturbed matter equation leads to 

- 2 (-iA'f{ + BiV'ifo) + f'{) + f' Q (-iwe^d + 8A[ - B[) + 2f ± (V"(J) - u?e 2A ) = (21) 

The gravitational equations are not independent : solving the second and the third for B\ , C\ and recon- 
structing C[ , B[ leads to an identity for the first equation. Knowing the expressions for B\ , C\ , B[ , C[ and 
further defining 

F:=h-^, (22) 
leads to the following master equation for the perturbation F(y) : 

F" - AA'F' + ( 4 t°ry 0) - 2 41- + + - V'ifo)) F = 0. (23) 



3A'(y) 9 A' 2 3 
Using the following change of function and change of variable: 

F -> ei A F, dy = e~ A dr, (24) 
the equation (1231) can be set in the form of a Schrodinger equation : 

2fjl_ 4^mi) lip _ 19 ? 2i - 
P 9A' 2 3A' 4 6 J0 J h K ' 

where we defined / (r) = f (y(r)), A(r) = A(y(r)), F(r) = F(y(r)). 

Note that the zero mode associated with the translational invariance of the background fields leads to a 
trivial F. Indeed the values A\ — A', fi — f imply F = 0. 



B. Fundamental solution 



For the range of the parameters considered, the fundamental solution is stable. Indeed, the potential is 
divergent close to the origin: 



A^A oa , /«/ ir =*V P «~, (26) 



r 



and the asymptotic behaviour is 



7T 



4V"(%) +15 



A«logr, /«- ^V p * ^ . (27) 

The profile of the corresponding potential is presented on Fig. [7] by the solid black line. The last term in 
the potential V p plays the role of the mass term, so it should be positive, leading to a positive fall-off of the 
eigenfunction. As a consequence, the potential is positive definite and there cannot be imaginary modes. 
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FIG. 7. Solid: potential of the perturbation master function for the solutions with one, two, three nodes. Dashed: 
the first normalisable eigenmode of the two node solution. 



The behaviour of the potential close to the origin ([21)1 and at infinity (|27p are the same as in the case 
of the fundamental solution, but the potential becomes negative in the intermediate region. The potential 
corresponding to the N = 2 and N = 3 solutions are represented on Fig. [7] by the red and blue lines 
respectively. The occurence of valleys where the potential is negative leads to bound states which possibly 
reveals the existence of instable modes if cj 2 < 0. Indeed, we find an unstable mode for the generic solutions 
with one mode. For k = 1.0 we find w 2 ~ —0.695. The profile of this eigen-mode (the normalisation is 
arbitrary) is supplemented on figure [7] by the dashed, red line. We further checked that the next bound state 
is indeed stable; this suggests that the brane characterized by N nodes of the scalar field has indeed N — 1 
unstable modes. This speculation was, however, not tested further. 



The purpose of this note is to exhibit new families of braneworld associated to five-dimensional grav- 
ity supplemented by a scalar sector admitting kink solutions. The family of solutions is organized like a 
tower of excited states of a Schrodinger equation (the bulk cosmological constant playing a role of spec- 
tral parameter, taking discrete values) although the underlying equations are highly non-linear. The types 
of braneworlds obtained in the paper have, to our knowledge, not been constructed yet and present some 
potentially interesting features for the localization of both species of chiral fermions. 

Often, when non linear equations admit various branches of solutions, the higher excitation are less stable 
than the main solution. It seems that this feature is obeyed in the present case : using a particular channel of 
the perturbed equations we checked that a typical solution with no nodes is stable and that the corresponding 
solution with one node is not. However, it could be expected that these new types of solutions is relevant 
in the mechanism of brane formations. Assuming brane formation in a dynamical process, possible unstable 
higher energy might be reached. The latter should then decay to the fundamental low energy configuration. 
These possibilities need further investigations. 



C. 



Excited solutions 



VI. CONCLUSIONS 
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